Astronomy & Astrophysics manuscript no. paper 


2nd February 2008 


(DOI: will be inserted by hand later) 





Spectroscopic diagnostic for the mineralogy of large dust grains 

M. Min, C. Dominik, and L. B. F. M. Waters 

Astronomical institute Anton Pannekoek, University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 
e-mail: mminOscience . uva . nl 

Submitted to A&A Letters 

Abstract We examine the thermal infrared spectra of large dust grains of different chemical composition and mineralogy. 
Strong resonances in the optical properties result in detectable spectral structure even when the grain is much larger than the 
wavelength at which it radiates. We apply this to the thermal infrared spectra of compact amorphous and crystalline silicates. 
The weak resonances of amorphous silicates at 9.7 and 1 8 pm virtually disappear for grains larger than about 10 pm. In contrast, 
the strong resonances of crystalline silicates produce emission dips in the infrared spectra of large grains; these emission dips 
are shifted in wavelength compared to the emission peaks commonly seen in small crystalline silicate grains. We discuss the 
effect of a fluffy or compact grain structure on the infrared emission spectra of large grains, and apply our theory to the dust 
shell surrounding Vega. 



1. Introduction 

Infrared spectroscopy is an invaluable tool for the study of the 
structure and composition of interstellar and circumstellar dust. 
In the infrared, functional groups in a dust grain lead to spec- 
tral features in the absorption cross section of a grain. These 
features can be detected in absorption against a strong infrared 
background source, or as emission features if the grains are 
warm enough to emit thermally in this wavelength region. The 
precise shapes and positions of such features contain informa- 
tion about grain size, shape and detailed chemical composition 
jBohren & Huffmanlfl983l) . 

It is generally assumed that dust grains need to be small 
in order to produce spectral structure. When grains are large, 
each individual grain becomes optically thick at infrared wave- 
length, and the features in the absorption cross section are ex- 
pected to weaken with increasing grain size and eventually 
to disappear completely. In most environments in interstellar 
space, small dust grains dominate the grain surface available 
and therefore also dominate the total absorption cross section. 
There are, however, environments where small grains are heav- 
ily depleted. In particular in gas-poor circumstellar disks, small 
dust grains are removed by radiation pressure and larger grains 
start to dominate the interaction with radiation. For example, 
the zodiacal dust in the solar system is dominated by large 
grains since grains smaller than about 1 micron are quickly lost 
tBurns et all Il979l) . This effect becomes even stronger in de- 
bris disks around A-type main sequence stars like Vega, where 
the increased luminosity shifts the typical blowout size to about 
10 micron jArtvmo wicz, 1988). Even stronger effects may ex- 
ist in disks around (post)-AGB stars, where already grains of 
mm size will be lost unless the disk is optically thick and gas- 
rich ( Dominik et all I2003I) . Another application could be the 



emission spectra of asteroids of which the regolith is depleted 
from small dust grains and references therein). 

We therefore examine in this paper the optical properties 
of large dust grains. In section[2]we show that the strong res- 
onances observed in crystalline materials lead to observable 
structure in the infrared spectra of large dust grains. In sec- 
tion we apply these results to the dust shell around Vega. In 
sections|4]and|5]we discuss and summarize our results. 

2. Infrared emission spectra 

In order to calculate the thermal emission spectrum of a dust 
shell we need to know the emission cross sections of the 
dust grains. With grain size, shape and chemical composition 
known, the emission cross section per unit mass as a function of 
wavelength can be calculated from the wavelength dependent 
bulk refractive index. In grains much smaller than the wave- 
length of radiation, a resonance i n the bulk refractive inde x 
causes a so-called surface mode (Bo hren & Huffman! IT983). 
These surface modes are strong enhancements of the absorp- 
tion cross section and are detected as emission features in the 
thermal emission spectra of these grains. For larger grains, 
the surface modes disappear and it is generally assumed that 
the thermal emission from large particles will be a smooth 
blackbody without a characteristic signature of the mineralogy. 
However, this assumption only holds if the refractive index is 
close to unity over the entire wavelength range. For the reso- 
nances in crystalline silicates, both the real and imaginary part 
of the refractive index can reach values up to ~ 9. For large 
particles with these values of the refractive index, most radi- 
ation will be reflected off the surface of the particle and not 
penetrate the grain. This leads to a significant reduction of the 
grain absorption efficiency at these wavelengths. Therefore, we 
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Figure 1. The emission cross sections per unit mass, k, for various sizes of amorphous olivine grains (left panels) and crystalline 
forsterite grains (right panels) as functions of the wavelen gth. The emissivit i es we re calculated using a uniform distribution of 
spheroids (UDS). The refractive index data are taken from Servoin & Piriou ( 1973). Note the different y-scales for the different 
grain sizes. 



also expect to find a corresponding decrease of the emission ef- 
ficiency, at wavelengths where there is a resonance in the bulk 
refractive index. This is illustrated in Fig.[Qwhere the emission 
cross sections per unit mass is plotted for grains composed of 
amorphous olivine and crystalline forsterite for different grain 
sizes. Moving from small to large grains, the spectral struc- 
ture in the spectrum of the amorphous olivine virtually disap- 
pears. The resonances in the spectrum of the forsterite grains 
change from emission enhancements to emission dips. These 
dips will be detectable in the thermal emission spectra from 
large compact particles. Generally the features at shorter wave- 
lengths and the strongest features are the first to disappear and 
change into emission dips. Most dips in the spectra of the large 
grains are slightly red-shifted compared to the corresponding 
emission enhancements in the small grain spectra. When the 
grain size is increased even more, the overall shape of the spec- 



trum will flatten but the spectral structure remains visible at 
about the same contrast level. 

In order to calculate the emission cross sections, a parti- 
cle shape has to be assumed. We used a uniform distribution 
of spheroids which is in very good agreement with labora- 
tory measurements of the mass absorption coefficients of small 

J F I 1 

crystalline silicate grains (Minetal., 2003b). The cross sec- 
tions w ere calculated using the method described in lMin et alJ 
(2003a) . For the large particles, the emission efficiencies are the 
same for all convex shaped particles. Calculations for homoge- 
neous spheres using Mie theory (|MieJ,|l908) will therefore give 
the same emission spectrum in the large particle limit. 

It is important to note that for large compact grains the 
emission cross section per unit mass is very low compared to 
small grains - this can be seen from the different scales used in 
figE Also, the emission dips created in large grains are weaker 
than the emission peaks due to small grains. In most spectra, 
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Figure 2. Mod el infrared spectra f or Vega. The parameters for the model are given in the text. The IRAS data are taken 
fromlWalker & Wolstencrofll dl988l) . the ISOPHOT data are from iHeinrichsen et all dl998h and the submm flux is taken from 
iHolland et al.Nl998l) . The labels with the curves in the inset are the abundances of forsterite for the different models. 



the emission by small grains, which is much more efficient, 
will be dominating. The spectral signature of the large com- 
pact dust grains should therefore be observable only in envi- 
ronments where small grains are heavily depleted. 

3. The spectrum of Vega 

An efficient mechanism to remove small grains from a popula- 
tion of grains is radiation pressure on grains in orbit around a 
star. Due to the efficient absorption, the radiative acceleration 
on small grains will easily exceed the gravitational force ex- 
erted by the stars. If the grains are orbiting in a gas-poor, opti- 
cally thin circumstellar shell, they will be removed on a Kepler 
timescale. Such conditions are given in the debris disks around 
main-sequence stars. 

To demonstrate the effects of large compact crystalline 
grains, we calculated the spectrum of Vega using a dust shell 
in a distance range of 80-120 AU from t he star. The s tellar 
spectrum is taken from a Kurucz model jKurucTi IT993I) with 
T e ff = 9500K and log g — 4. For the size distribution of the 
dust grains we used a powerlaw with sizes between 15 and 
100 fim and a powerlaw index of -3.5, typical for a collisional 
cascade (Tanaka et al., 1996). From the collisional model of 
iDominik & Decinl d2003l) and the calculated radiation pressure 
we estimated the fraction of small particles to be < 0.1% of 
the total mass. Calculations show that this mass fraction does 
not produce significant emission to be visible in the spectrum. 



The dust is assumed to be a mixture of 20% amorphous carbon 
and 80% olivine by mass. We calculated various models with 
different mass fraction of crystalline and amorphous olivine. 
The dust grains are assumed to be compact particles, the opti- 
cal properties are calculated using homogeneous spheres. The 
dust temperatures are computed self-consistently from the local 
radiation field. The total dust mass used in order to fit the pho- 
tometry points is ~ 1 /TO of a lunar ma ss which is in agreement 
with the estimate bv bent etal] (E000). 

The model spectrum is shown in Fig.|2] which runs reason- 
ably well through the observed photometry points. It should 
be stressed here that the presented spectrum is an example to 
demonstrate the effect of large grains. We do not claim to pro- 
vide the best fit to the available data points. The inset in Fig. [3] 
shows an enlargement of the 20 to 50 /im region. A promi- 
nent dip in the spectrum at 34 jjm is caused by large com- 
pact forsterite grains. Even with only 5% of the mass present 
as compact forsterite grains, the dip would most likely be de- 
tectable by SIRTF. 

4. Discussion 

4.1. Grain composition 

In section[3]we have shown that a fraction of about 5% in com- 
pact forsterite grains would produce a detectable feature in the 
spectrum of Vega. It is important to note that this mass frac- 
tion should be present as seperate homogeneous grains. If the 
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forsterite material is present as small inclusions inside an amor- 
phous matrix, a significantly larger fraction of forsterite would 
be required to produce a spectral signature. 

4.2. Grain structure 

The structure of the dust grains is an important parameter that 
contains information about its formation and processing his- 
tory. As a first rough classification we can distinguish between 
fluffy and compact particles. The structure of a grain can be 
traced back to its formation process. Dust grains that are grown 
by coagulation of small (submic ron sized) gra ins will show a 
very fluffy open structure dBlum et all l2000l) . similar to the 
structure observed in interplanetary dust particles (IDP). Such 
grains are also expected to result from the destruction of prim- 
itive planetary bodies such as small cometary nuclei. Compact 
(non-fluffy) large grains will result from the destruction of 
larger parent bodies which have been molten and possibly dif- 
ferentiated during the formation process, like large asteroids. 

Measurements of the infrared t ransmission spectrum of a 
large (a * 20 /mi) fluffy IDP by iMolster et all J2003h still 
shows the resonances characteristic for submicron grains. This 
indicates that the infrared spectrum of fluffy particles will still 
be dominated by the spectral signature of the constituents of 
the aggregate. In other words, we expect to see the small parti- 
cle crystalline silicate emission features in the emission spectra 
of these dust grains. For large compact particles we expect to 
see a spectral signature as shown for 100 /mi sized particles in 

Fig.m 

5. Conclusions 

We have developed a new diagnostic tool for determining the 
mineralogy of large compact dust grains. The thermal emission 
spectra of large compact crystalline silicate dust grains exhibit 
characteristic features that can be detected with high signal to 
noise spectroscopy. These features are emission dips at posi- 
tions that can be slightly red shifted from the corresponding 
emission peaks observed for small grains. As an example, a 
model calculation of the infrared emission from the dust shell 
around Vega shows that these features should be detectable 
if the dust grains are compact and contain at least 5% crys- 
talline forsterite in a seperate grain component. Since there is a 
clear spectroscopic difference between fluffy and compact dust 
grains, the spectral signature of crystalline silicates can be used 
to obtain information on the structure, mineralogy and forma- 
tion history of the dust grains. Other applications of this diag- 
nostic can be the zodiacal dust or asteroidal regolith. 

Our model calculations show that if large compact crys- 
talline silicates are present even with modest abundances, this 
should be detectable as a depression near 34 /mi. If the grains 
are fluffy, an emission band should be present at a slightly 
shorter wavelength. Such spectral structure may be measured 
using the IRS spectrograph on board of SIRTF. 
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